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Abstract 

We present new improved parton distributions for the photon. We fit all 
available data on the photon structure function, i<2^(x,Q^), with > 3 
GeV^, in order to determine the quark distributions. We also pay particular 
attention to the gluon distribution in the photon, g'^{x,Q'^), which has been 
poorly constrained in earlier analyses which only include structure function 
data. We use large pT jet production in 77 collisions from TRISTAN to 
constrain g'^ . We also see what information can be gleaned from ^p collisions 
at HERA on g'^ and on the quark distributions at large x, where there are no 
structure function data. We review future prospects of elucidating the parton 
distributions of the photon. 
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I. INTRODUCTION 



Some time ago, we presented some parton distribution functions (pdfs) for the photon 
These distributions, which we will refer to here as the GS distributions, were the first 
available in both leading order(LO) and next-to-leading order(NLO) QCD. Since then they 
have been used as input to perturbative QCD (PQCD) comparisons to a variety of types of 
data with a fair degree of success. Applications in 77 collisions at various e"'"e~ colliders are 
not restricted to the analysis of F2{x, Q"^) data 0; they also include the inclusive production 
of jets and single hadrons at large px Applications at the ep collider HERA include 
the production of jets and dijets PHTSl and single hadrons at large pt- Up to now, 

no contradiction with any data has been found. These pdfs have also been used in estimates 
for rarer processes such as the photoproduction of high-mass di-lepton pairs and large 
Pt prompt photons |T^JT8[] . 

At that time, and up until fairly recently, the only experimental information on the pdfs 
of the photon was obtained from studies of the structure function of the photon, F2{x, Q^), 
in two photon collisions at e~^e~ colliders. In LO QCD, F2{x,Q'^) is related to the quark 
and anti-quark distributions in the photon, ql{x,Q'^) and q]{x,Q'^) respectively, by, 

Nf 



— T 

i=l 

■)2\ 



Fnx,Q')=T.^e^{<in^,Q') + qnx,Q')) (1.1) 



and thus the quark distributions are fairly directly determined from the data on 
Theoretically, this is a very interesting area, because, at large x and asymptotically large 



Q^, the pdfs of the photon and hence F2 (x, Q^) are predicted from PQCD in both LO [19 



and NLO [^^. The so-called 'anomalous' piece is calculable as a function of x and 
and it is the dependence (specifically oc (a^ 

(g2))-^) that ensures that it dominates 
at asymptotic values of Q^. However, it is singular as x — * and to regularise these 
singularities, some non-perturbative input to the Altarelli-Parisi(AP) equations is required, 
at a reference = Q^. The anomalous {asiQ"^)) ^ behaviour of the pdfs arises because of 



the direct 'J ^ qq coupling, which gives inhomogeneous terms in the AP equations [22|: this 
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is the special feature of the photon. However, in the range experimentally accessible at 
present and in the foreseeable future, the effect of this input at a = significant, and 
so the analysis of the photon structure function is similar to that of a hadronic structure 
function. Different groups make different choices of input scale different forms of input 
pdfs, and all fit parameters to Q"^) data at various Q^. The different possibilities have 

been recently reviewed in ref |^ . 



The gluon distribution of the photon, g^{x, Q"^), are not well constrained by these struc- 
ture functions analyses [0, essentially because the coupling of the AP equations for the 
gluon and singlet quark sectors is weak, and so the output (evolved) quark pdfs and thus 
the evolved values of F2{x,Q'^) (for Q"^ > Ql) are rather independent of the input gluon 
distribution at = Qq, g^{x,Ql). Furthermore, and even more importantly, unlike the 
case for hadrons, g'^ is not constrained by a momentum sum rule [0,0]- The result of all 
this is that the available parametrisations of the photon distributions agree reasonably well 
in their quark distributions but can have considerably different gluon distributions p3|-p5|. 
The input gluon distributions in the evolution equations, while not completely arbitrary, are 
currently just theoretically motivated guesses. 

Since our distributions were published, there has been much experimental activity in this 
and related areas. Not only has the amount of data on F2 been greatly increased, but data 
on jet production in both 77 collisions at e~^e~ colliders and in jp collisions at HERA have 
appeared. These provide information on both quark and gluon distributions in the photon 



via resolved photon processes ||2^, as we shall see. 

In view of this and also in view of the fact that there is expected to be an explosion of 
data in all of these areas in the next few years, it seems like an opportune moment to provide 
improved versions of our pdfs. The new data will come from higher luminosity running at 
HERA and also the commissioning of LEP2, which is anticipated to provide several hundreds 
of inverse pb of integrated luminosity in e~^e~ collisions over its lifetime. The improvements 
we anticipate making follows: 

(1) Extending the kinematic range of validity in from > 5.3 GeV^ down to > 3 



(2) Fitting to all available data on to obtain a better determination of the 
quark distributions. 

(3) Constraining the input gluon distribution by fitting the TRISTAN data on jet pro- 
duction in 77 collisions. We will also see what can be learned from HERA data on jet 
production. 

The plan of this paper is as follows. In subsect. II. A we discuss general considerations 
about the pdfs of the photon. In subsect. II. B, we set up the parametrisation of the LO pdfs 
and in subsect. II. C we discuss our prescription for going over to the NLO parametrisation. 
In both cases we pay particular attention to the large x behaviour. In sect. Ill we present our 
fits to the data, considering separately the data on F2 (subsect. III. A), jets in 77 collisions 
(subsect. III.B), and jets in ■yp collisions (subsect. III.C). In sect. IV we discuss the 
properties of the new distributions and review future prospects for improving our knowledge 
of the photon pdfs. In sect. V we present conclusions. 



II. NEW PHOTON DISTRIBUTIONS 



A. General Considerations 



To obtain pdfs for the photon, we must start with input photon pdfs at some reference 
scale Ql, and choosing an appropriate value of Ql provides the first problem. Making a 
vector meson dominance (VMD) ansatz seems an obvious starting point as VMD provides 
a connection between the photon and p meson pdfs. If we then use SU(6) to relate the p 
pdfs to those of the pion, we can then use experimental constraints on the pion pdfs from 



Drell-Yan lepton pair and prompt photon production data [^. In such a VMD ansatz, the 
gluon is fixed by a momentum sum rule via a mesonic momentum sum rule. However, a 
problem arises here: if we choose a plausible VMD scale, — 1 GeV^, a pure VMD input 
is known to be insufficient to fit the data at higher [28- p0|] . To circumvent this problem 
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two separate approaches have been adopted. The first is to maintain the VMD approach 
but to start the AP evolution at a very low scale significantly below Qq — 1 GeV pT|-p^. 
The other is to take an input scale significantly above Qo — 1 GeV and fit the quark pdfs to 
the F2 data here, which essentially means supplementing the VMD input with a point-like 
component, seemingly naturally provided by the Born-Box diagram [^^^|-|^. Unfortunately 
there is no correspondingly natural choice for the gluon density and a guess must be made. 
This is the approach that we will adopt here, as we did in ref , though we will be reducing 
from the value chosen there, Ql = 5.3 GeV^, to Qo = 3 GeV^. In ref [IH we advocated 
choosing the gluon content to ensure that the ratio of momentum carried by the quarks and 
antiquarks relative to that carried by the gluons, R2{Q^), should lie between 1 and 3. The 
motivation for these limits is as follows. At low we expect VMD to be good, the photon 
to be hadronic, and i?2 — 1 as for a hadron: at high Q^, we expect the asymptotic result 
of PQCD |3^, that i?2 = 99/32 ~ 3 to hold, and so we expect a steady increase from 1 to 
3 as increases. It is interesting that at the time we noted that the LAC distributions 
have such a large gluon component that they are in danger of violating this bound. Since 
then LAC3 has been rejected by the data, because the large gluon component at large x 
overestimates the jet cross section in both 77 and 7p collisions and LACl seems to have too 
many gluons at small x compared to the HI data. We will expand on this later. 



B. Input Distributions: LO case 

To obtain regularised distributions in LO we separate the photon structure function into 
a hadronic and a point-like part p8[ 

F^ix, Ql) = Fi^x, Ql) + Ql) (2.1) 

The hadronic part, which we assume can be approximated by VMD, is important at small 
values of x (< 0.2) while the point-like part, which we base on the Quark Parton Model 
(QPM) formula, models the medium to large x part. 
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The hadronic part of the input was chosen according to standard VMD ideas. We assume 
that the hadronic photon can be represented by the p° meson 

FlvMnix,Q') = ^E^'^^fi^^Q') (2-2) 

J p i 

where fp/4:7i ^ 2.2 , and use SU(6) and isospin invariance to relate the p° distributions to 
those of the charged pions which are constrained by data [ p7[] : 



Qf{x,Ql) = qf{x) = liqfix) + qf {x)) (2.3) 

The above assumptions are supported by the fact that the VMD estimates of the structure 
function F2{x,Q'^) base agree with the data at Q"^ = 1 GeV^ HSJSO]. We also make the 



plausible assumption that the VMD breakdown into singlet, non-singlet and gluon sectors 
is also reliable. In terms of valence v^{x) and sea, C"{x), these are given by 



Ana 



^'vMDix, Qt) = ^^^i^i^v- + 6C) (2.4) 



Ana 



p 



Jp 

where we have used Nf = 3 flavours since the heavy quark contributions are expected to 
be very small at this Q^. The constant k, where 1 < /t < 2, was introduced in ref. ||3^ 



to represent the uncertainty over the inclusion of higher mass vector mesons. We take the 
following forms [jl[: 

Xf '^(x) = Ay^{l — x) 

xC{x) = B{l-x)^ (2.5) 
xG^{x) = C{l-x)^ 

The conservation of baryon number fixes A at 0.75, fixing the momentum fraction of the 
hadronic part carried by the valence sector at 40%. The momentum fractions carried by the 
sea and gluon sectors are fixed by B and C which we use, along with k as free parameters 
in our fits to the Fn data. 



For the point-like part we use the lowest order Bethe-Heitler (B-H) form pO| , p 



vr 



2 



(8x(l - x) - 1 - - x) 



for > where 



(2.6) 



f3^ = 1 



:i-a;)g2 

and irii are the light quark masses which are taken as fairly closely constrained parameters. 
Now q]{x, Q2) = for /^^ < 0, i.e. for 



g2 + Ami 



l>x>x, = T^f^2 (2-7) 



and typically ~ 0.9 for light quarks and Qq = 3 GeV^. Also the VMD piece vanishes as 
X — s> 1 as a power of (1 — x). Hence our quark distributions will be greatly suppressed at 
large x, a point we will return to later. The parametric form for the distributions are thus 

ql{x,Ql) = q^'^ix.m^.m,) ^ qf^^'ix.K.B.C) (2.8) 

for the singlet and non-singlet sectors and for the gluon we take the form: 

9\x. Ql) = 9{x, K, C) + |-AP,,(x) * S(x, mO (2.9) 

Po 

and the convolution * is defined by 

A{x) * B{x) = —A{x/y)B{y) (2.10) 
Jx y 

and A is a parameter which enables us to adjust the component of the gluon input estimated 
from bremsstrahlung off the singlet quarks. This will increase the fraction of the photon 
momentum carried by the gluons, and compensate for the increase in the momentum carried 
by the quarks due to inclusion of the point-like contribution. In our earlier work where 
we took Ql = 5.3 GeV^, we presented two distributions, GS(I) and GS(II), corresponding to 
A = 1, respectively: these numbers cannot be compared directly with those in the present 
work, different value of Ql is used. 



C. Input Distributions: NLO case 



The input distributions in NLO were obtained from the LO ones by demanding that the 
structure function, F2 (x, Qo), be the same at Qo = 3 GeV^ in LO and NLO i.e. 

F',''''^{x,Ql)=F',^'^''Ml) (2.11) 

in the same way as in ref. [Q]. For the gluons we simply take 

but choosing the NLO quark distributions in order to satisfy the above condition is a non- 
trivial procedure which highlights the well-known fact that pdfs are renormalisation scheme 
dependent and have no physical significance on their own. 

To proceed, we assume that the distributions in NLO can be obtained from the LO ones 
by the addition of a correction term 

g7'^^^(x,Qg) = g7'^''(x,Q^) + aM)<i\^.Ql) (2-13) 

These are substituted into the defining equation for F2 (x, Qq) which is 

^f:j{x, Q') = (^5(1 - x) + ^f!^B,ix)j nlsi^, Q') 

+ <e^> [5{l ~x) + ^^!^Bg{x)j Q2) (2.14) 



<e'> ^^BG{x)*g^x, Q') + 6,B,{x), 



where 



5^ = 3iV;<e^>^. (2.15) 
47r 

and the functions Bg{x), Bg{x) and B^{x) are the Wilson coefficient functions. as(Q^) is 
now the two loop couphng constant defined by 



, 9. 4:71 



/3oln(QVA2) 



^ Aln[ln(QVA2)] 



/3o^ln(QVA2) 



(2.16) 



Combining these equations at Q"^ = Ql with a few extra assumptions to get the singlet and 
non-singlet distributions enables us to extract the q^^^ and hence NLO quark distributions 
Of course, to carry this out we need to specify the renormalisation scheme as the Wilson 
coefficients and hence the pdfs are scheme dependent. It is only in combinations such as 



eq.(2.14) that they have physical significance. We work in the MS scheme: expressions for 
the Wilson coefficients in this scheme are given in ref . 



There is a subtlety peculiar to the photon in the MS scheme which is discussed in detail 



in ref. [^p3[ . It arises because is divergent and negative as x — 1 {B^ oc ln(l — x)). This 
means that if we put our LO input into eq.(2.14), we would get a negative -^^(x, Ql) at large 
X, though the problem goes away at very large Q^. The method outlined here eliminates this 
problem of negative values of F2{x,Ql) at large x due to the B^{x) term in the definition 
of the structure function. What we have done in imposing eq.(2.11) is essentially to add a 
term to the quark distributions to compensate for the B^ term. The resulting distributions 
are steeply rising with x and are quite different from the LO ones, as we shall see. 

III. FITS TO DATA 

A. Structure Function Data 

Here we fit all the data on with > 3 GeV^: we include data from PETRA g^-gSf, 
PEP H-P, TRISTAN [IHH, and LEPl ||,|6l. The data cover a wide range of and 



some have the charm subtracted out and some do not and so we must have a consistent 
policy over its treatment. For Ql < < Qj = 50 GeV^ we use Nf = 3 evolution and 
when we compare to data which includes charm we add to F2 the B-H estimate of charm 
as provided by eq.(2.6), taking = 1.5 GeV. For > Ql, all of the data have charm 
included. For these, we use a Nf = 4 evolution which has been started from = Q2 = 10 
GeV^ with inputs generated by our Nf = 3 evolution from Ql to Q2 supplemented by the 
addition of the B-H term at Ql- Throughout we keep fixed at 200 MeV. 
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In figs. (1-4) we show the fits to the data on the ^-distributions at various in both LO 
and NLO. To confirm the insensitivity of structure function analyses to the gluon content 
of the photon we show a sample of fits corresponding to different values of A. It can be seen 
that the curves are indistinguishable, except at small x. 

In fig, 5 we show the fit to the data on 

r\^{x,Q^)dx (3.1) 

JO. 3 

plotted against ln((5^) for LO and NLO. The two curves agree at Qq but it can be argued 
that at higher Q^, the NLO evolution gives a somewhat better agreement with the data. It 
should also be noted that our method of heavy flavour inclusion, discussed above, leads to 
no discontinuity in the curves when one passes the heavy flavour thresholds. 



B. Jets in 77 collisions 

Until relatively recently, studies of jets in 77 reactions were confined to model-dependent 
analyses of multiparticle distributions such as thrust. These were able to prove the existence 
| p7| , |55[ ] of the contribution from resolved photon processes but were not precise enough 



to distinguish between different photon pdfs, though some effort was made along these lines 
by the later analyses indicating some preference for the GS [0 pdfs. 

Since then data on the production of jets in appeared from the TRISTAN collaborations 
AMY |^0[ and TOPAZ [Q,^]. Here we will attempt to constrain the input gluon distribution 



by fitting to their data on single and two-jet production. These have been shown to be 
somewhat sensitive to g'^ in that they cannot be fitted with g"^ = and also have been 
used to rule out the somewhat extreme LACS parametrisation, which, with its large 
gluon component at large x, gives much too high a cross section. Some calculations of 
jet production in 77 collisions in NLO have been carried out P,|6|,|62| , which have mainly 
concentrated on the theoretical uncertainties, although some preliminary phenomenology 
has been done In view of the theoretical and experimental uncertainties we feel that 

at the moment it is only worth carrying out a LO comparison with the data. 
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In this, for the photon flux we use the Equivalent Photon Approximation(EPA) |63|, with 
parameters fixed according to the experimental conditions. We keep our quark distributions 
fixed at the values we obtained in subsect. IIIA, and fit the data by varying the gluon 
content achieved by varying the parameter A. We find good fits to the data with A in the 
range 0.90-1.0. We show our fits to the data in figs. 6 and 7. To illustrate our sensitivity to 
the gluon we also show sample fits with g"^ = and those corresponding to A = and 1 



C. Jets in ^yp Collisions 

At HERA, ep interactions in the untagged case are dominated by photoproduction with 
the spectrum of the (mainly) real photons given by the EPA [Q. The first data on jets ||64 



were superseded by the first jet cross section data which enabled the direct photon reactions 
to be separated from the resolved photon reactions and for the LACS ||36| pdfs to ruled 



out as they overestimated the jet cross section because of their large gluon distribution. 



Since then much better quality data on jet and di-jet cross sections |TT] , p^|67| , |68| have been 
published which it is hoped can be used to give information on the pdfs of the photon. 
These data have been confronted by NLO calculations i69|j70| j9[| and the fits are not good 
| T0| . What conclusions can be drawn from this is not clear at the moment as there are other 



aspects of the data which are difficult to understand, for example, the energy flow out of 



the jet for the resolved component [^]. Multiple parton interactions have been suggested 



as the reason but it is not clear how they can be incorporated into a NLO calculation. 
In a LO calculation, they give an effect in the right direction, increasing the cross section in 
the regions which are dominated by the resolved contribution, particularly at lower pt, and 
they also help with the energy flow problem. 

In view of this confusion, it seems premature to compare the results of a full NLO 
calculation with the data: we will simply confine ourselves to a few semi-quantitative remarks 
on what has been learned from the data that is relevant to pdfs. The demise of LACS shows 
that the we cannot have a large gluon distribution at large x. There are two other places 
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where something relevant can be learnt. 

The first is the fact that many pdfs overestimate the photoproduction cross section 
in the region of extreme negative rapidity. This is true for both single jet and 
dijet p.|TO[] cross sections. Here the cross section is supposed to be dominated by the direct 
contribution which depend only on the proton pdfs which are well-known from deep inelastic 
scattering from protons. In direct photon processes multiple parton interactions cannot be 
the culprit: in fact the energy flow is well understood here. It turns out that with the NLO 
pdfs used initially GRV and AFG, the resolved contributions make a significant 

contribution here, presumably because their quark pdfs are large at large x, and a parton 
with large x contributes like a direct photon, as it brings all of the photon energy into the 
hard subprocess. This problem goes away if the GS distributions are used both in the dijet 
case P,p!0| and in the single jet case |T^. This is presumably because they are relatively 



suppressed at large x, as we pointed out earlier, certainly in LO where the pdfs have physical 
significance, and this relative suppession is reflected in the NLO quark distributions p3[ . 



This difference in the quark distributions does not show up in structure function analyses: 
there are no data at large x and are unlikely to be in the future as large x corresponds to 
W, the energy of the 77 — > hadrons, process being very small. 

The other place is the gluon distribution of the photon which has been "extracted" by 
HI ||6^, They used the following procedure: they constructed an x^ distribution of the dijet 
cross section, subtracted the contribution from direct photons to obtain the resolved cross 
section, and then subtracted the quark contribution from that, leaving the gluon distribution 
of the photon, which we plot against ours in fig. 8. There are many uncertainties associated 
with this procedure: it is only defined in LO , scale dependence is obviously a problem, 
as are multiple interactions, which are larger at small x^. Also the quark contribution 
subtracted was that corresponding to GRV: this is probably irrelevant in view of the other 
uncertainties catalogued above. However, the result is perhaps indicative that the data 
prefer more moderate gluon distributions such as ours and GRV. It is the extreme gluons 
such as LACl and LACS that have problems. 
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IV. DISCUSSION 



The new GS photon distributions are many respects similar to the GSl set of ref |l|]. As 
we have discussed, the main improvements are; a lower starting scale for the evolution and 
more importantly, we have fitted to all available data on the photon structure function as well 
as attempted to constrain the gluon input using single and two jet data from TRISTAN. The 
HERA data on single and two jet production has also provided some indication that smaller 
quark distributions such as ours are preferred at large- a; and not the large distributions 
given by most of the other parametrisations on the market. We have attempted here to use 
information from all possible sources, both recent and older, experimental and theoretical, 
in order to constrain the parton distributions of the photon. For example the ratio R2{Q^), 
discussed in the introduction, which is the ratio of the momenta carried by quarks to gluons 
starts at around 1 and approaches 3 as increases. We regard this as a kind of theoretical 
constraint on the gluon distributions. 

There are some areas where some improvements could be made in our parametrisations. 
Two such areas are the flavour decomposition of the structure function, and continuity 
in of the quark distributions. In our treatment, the four flavour evolution is started 
at = 10 GeV^ but not used except above = 50 GeV^. At this point the charm 
distribution is taken as equal to that of the w-quark, which undoubtedly overestimates it by 
a considerable amount, but correspondingly, the w-quark distribution is reduced at this point, 
rendering it discontinuous there. The sum of the quark distributions, and hence F2 {x, Q"^) 
on the other hand has no discontinuity at this point. We therefore recommend that single 
flavour distributions, such as required for charm production in a massless framework not be 
calculated using our distributions. 

As mentioned above, although there are no data on Q"^) at large-x, above x ~ 0.8, 

the recent HERA data on jet production appears to prefer smaller quark distributions in this 
region. At small-x, below x ~ 0.05, on the other hand we have no experimental information 
and hence cannot constrain the distributions here, except via the theoretically motivated 
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models. It is hoped that at in the near future this situation may be remedied. 

Our knowledge of the photon distributions has increased very significantly over the past 
few years, thanks to new experimental data from TRISTAN, LEP and HERA. The future 
outlook is even more optimistic, as HERA is expected to run at higher luminousity and 
LEP2 will soon start taking data. We anticipate that processes such as prompt photon | |73[| , 
heavy flavour and lepton pair production will be measured with good statistics at these 
machines, allowing tighter constraints on the quark and gluon distributions of the photon. 

V. CONCLUSIONS 

We have attempted to update the GS photon distribution functions by, among other 
things, constraining the input gluon distributions by fitting to TRISTAN jet data. We 
find that although the data is useful as an indicator of the existence of g"' and to rule out 
extreme distributions, it is not good enough to fix the distribution better than current VMD 
estimates, or other reasonable fits. 
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Figure Captions 

Fig. 1: Fits to PETRA data on tlie plioton stucture function. 
Fig. 2: Fits to PEP data on the photon structure function. 

Fig. 3: Fits to TRISTAN data on the photon structure function. At = 390 GeV^ in 
the second graph the solid curve is for the choice A = 2.0 (see text). 

Fig. 4: Fits to LEPl data on the photon structure function. At = 12 GeV^, in the 
second graph, the dot-dashed curve is for the choice A = 0. 

Fig. 5: dependence of the structure function integrated over the region 0.3 < a; < 0.8. 

Fig. 6: Fits to AMY data on 77 —>■ jets for two choices of A and for g'^ = 0; (a) fits to 
the pt distribution of the single-jet data (b) fits to the rapidity distribution of the 
single-jet data and (c) fits to the pt distribution of the two-jet data. 

Fig. 7; Fits to TOPAZ (a) single-jet and (b) two-jet data. 

Fig. 8; Comparison of the GS(HO), GS(LO) and GS(1) and GS(ll) LO gluon distribu- 



tions with the gluon distribution extracted by the HI Collaboration ||68 |. 
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